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Symbols 

a - orbit semimajor axis 
E - orbit eccentric anomaly 
e - orbit eccentricity 
F - thrust magnitude 
I - orbit inclination 

J 20 - coefficient in earth gravitational potential ex- 
pansion related to polar oblateness, 1.082 x 10“ 3 
J 22 - coefficient in earth gravitational potential ex- 
pansion related to equatorial oblateness, 

-1.8 x 10' 6 

M - orbit mean anomaly 
m - satellite mass 
n - satellite mean motion 
p - orbit semilatus rectum 

R,S,W - perturbing accelerations with respect to orbital 
reference acting on satellite 
R - mean equatorial radius of earth 
r - magnitude of satellite radius vector 
r s - earth-synchronous radius 
t - time 

V c - circular velocity, V c = na 

X,Y,Z - perturbing accelerations with respect to rectan- 
gular coordinates 

x,y,z - rectangular coordinates of satellite with respect 
to earth 

e - related to epoch by e = M - nt + o> 

6 - orbit true anomaly 
M - earth gravitational constant 
H - longitude of ascending node of satellite orbit 
ai - longitude of perigee of satellite orbit, ?/> = 0. + oj 
co - argument of perigee of satellite orbit 

Abstract 

Satellite orbit adjustment and station keeping can be 
accomplished by electric propulsion systems. The sys- 
tems considered in this study are those which are limited 
to thrust levels of a few millipounds. These thrust levels 
result in relatively slow variation of the state variables, 
and therefore the equations of motion lend themselves 
quite well to perturbation techniques of solution. Under 
these conditions the use of electric propulsion for orbit 
adjustment and station keeping and attitude control of 
typical earth-synchronous and sun-synchronous satellites 
has been studied. For a sun-synchronous meteorological 
satellite it was found that a five millipound resistojet sys- 
tem under continuous thrust application could correct an 
orbit precession error of 0.029 degrees per day by a 
spiralling operation lasting 11.3 days. It was found that 
the method of inclination adjustment and constantly bal- 
ancing the precession moment were less efficient than the 
spiralling approach. For a 1400 pound earth-synchronous 
communications satellite a five pound hydrazine system 


was assumed for removal of gross injection errors. A 
17.4 millipound resistojet was employed for vernier orbit 
adjustment of 22 feet per second which required 26 days. 

A self optimizing adaptive technique was assumed for at- 
titude control using a multijet thruster module system. 

Introduction 

Recently, the development of electric propulsion has 
been followed more and more attentively by designers of 
our future spacecraft. As a result, much work has gone 
into analysis of the application of electric propulsion to 
specific missions. In the first section of this paper the 
equations of motion are briefly reviewed for a satellite 
under the influence of small perturbative accelerations. 

The equations presented are those of Gauss, which are de- 
rived in Chapter XI of Reference 1. They were first used 
by Gauss in his studies of the perturbations of Jupiter 
caused by Pallas. These equations are then specialized 
for analysis of sun-synchronous and earth-synchronous 
missions. The purpose of this paper is to study these mis- 
sions for which electric propulsion is especially attractive. 

Among the better known of the sun-synchronous satel- 
lites are the Tiros series of meteorological satellites. In 
fulfilling its mission, this type of satellite is designed to 
provide meteorologists and other scientists with observa- 
tions of cloud cover and cloud-top temperatures, solar 
proton density measurements and heat balance measure- 
ments. The satellite being studied is required to be stabi- 
lized about the roll and yaw axes, and motion about the 
pitch axis is required to be controlled to provide an earth- 
oriented platform. A magnetic torquing control system, 
References 2 and 3, is well suited for the sun-synchronous 
satellite. However, the control system is limited because 
the torquing coil which generates the magnetic field within 
the satellite is limited in size. As with all satellite orbits, 
injection errors are inherent, as are misalignment errors 
of thrusters for orbit correction. Any disturbance torques 
resulting from orbit correction must be within the control 
system limitations. For this reason low thrust, and there- 
fore electric propulsion, is attractive. The relatively long 
time required for correction using electric propulsion is 
acceptable because the errors resulting from imperfect 
injection have slow accumlative effects. 

There have been recent significant steps taken in the 
development of spacecraft antennas. Progress made in 
antenna technology is one of the major factors which pace 
the advancement of on-board radio communications and 
radar in orbiting satellites. Many of these satellites are 
flown in an earth-synchronous orbit. As progress is made 
in antenna technology, progress must also be made in con- 
trol system technology for controlling the satellites. Typ- 
ical of the requirements are a fine pointing mode, to within 
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0.1 degrees of any point on the earth's surface, and a 
slewing capability, 17.5 degrees, to cover the earth's 
disc. Self optimizing adaptive techniques, utilizing elec- 
tric propulsion, are being studied to provide this control. 
For obvious reasons, it is advantageous to use the same 
propulsion system for as many functions as practical. 
Initial attitude acquisition, from relatively low body rates, 
and vernier orbit correction are two of these functions. 
Another function is that of station keeping or offsetting 
the perturbation caused by the non-circularity of the 
earth's equator, References 4 and 5. 

Equations of Motion 

General Equations 

The derivation and use of the equations of motion for 
one body under the gravitational influence of another pri- 
marily, but with "perturbing" accelerations from other 
sources, have been well documented in the literature. 

A brief history is given at the end of Chapter XI of Ref- 
erence 1. The purpose here is to review and extract 
those equations which are of direct interest in a study of 
orbit adjustment and station keeping for earth-synchron 
synchronous and sun-synchronous satellites. For details 
of the development the reader should refer to Chapter XI 
of Reference 1. 

For a satellite orbiting the earth, the equations of 
motion can be written as: 
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Specialized Equations 

These equations have been called Gauss' equations as he 
was the first to derive and use them. The three mutually 
perpendicular components of perturbing acceleration are 
defined as: R is the component in the direction of the 
radius vector (positive in the direction of increasing 
radius vector), S is the component perpendicular to R 
in the orbital plane (positive in the direction of increasing 
longitude in the orbit), and W is the component perpen- 
dicular to the orbital plane (positive in the direction in 
which the orbital motion appears clockwise). 

The general perturbation equations of the preced- 
ing paragraph can be specialized for analysis of sun- 
synchronous and earth-synchronous satellites. The orbits 
being analyzed are nearly circular and therefore the ef- 
fect of eccentricity is neglected, hi addition the differen- 
tials will be replaced by variations as the changes are 
relatively small. Equation 2-a can then be rewritten as: 


where the influence of the gravitation of a spherical earth 
is included on the left side. The x , y , and z are rec- 
tangular coordinates of the satellite with respect to an 
earth fixed coordinate frame. The X, Y , and Z are 
"Perturbing" accelerations due to the earth's oblateness, 
possibly the gravitational attractions of the sun and moon, 
and that due to thrust. After some rather lengthy but 
elegant mathematical manipulation by variation of con- 
stants there results: 
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where a Q is the semimajor axis from which the variation 
Sa is taken and the thrust F is applied in the direction 
of S. V is the equivalent circular velocity, V c = na. 

Under the same conditions specified above, small var- 
iations in inclination and motion of the node (precession 
rate) reduce from Equations 2-c and 2-d to 
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where in this case the thrust F is applied in the direction 
of W. Now if in Equation 4-b thrust is applied in the pos- 
itive direction from 6 = - a> to 9 = - w + n and is then 
reversed in direction from 8 = - a> + n to 8 = - w + 2-n 
an average rate of change of precession rate results 
giving; 


sin I 5 
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Similarly, when the change in direction of thrust is shifted 
shifted by 77 / 2 an average inclination change per orbit, 
from equation 4 -a, results as 


51 = - — ) St 

77 V 

c 

Equations 3, 5, and 6 will be of particular value in the 
following sections. 


coordinate system is assumed with origin at the center of 
mass and axes oriented as follows: 

X-axis or yaw axis - normal to satellite spin axis and 
surface desired to face the earth. 

Y-axis or roll axis - normal to satellite spin axis and 
pointing along velocity vector when satellite has de- 
sired orientation. 

Z-axis or pitch axis - lies along satellite spin axis and 
completes right hand coordinate system. 

Representative moments of inertia for this configuration 
are: 

I x = 633 inch-pound-seconds squared 

I = 570 inch-pound-seconds squared 

I = 562 inch-pound-seconds squared 


Sun-Synchronous Satellite 


Satellite Injection and Orientation 


Satellite Description 

The configuration of the sun-synchronous meteoro- 
logical satellite being considered is shown in Figure 1. 
The orientation is that desired when the satellite is in its 
final operating configuration. The main features are the 
equipment module (main body) and the deployed three- 
panel solar array. The central equipment module or 
main body is a rectangular prism. The base of the prism 
is approximately 40 by 40 inches. This module houses 
the data gathering and associated support subsystems. 
Also housed within this module is a momentum flywheel 
which provides gyroscopic stabilization of the satellite. 
The satellite angular momentum vector is assumed to 
have a magnitude of 200 inch-pound-seconds. 

A weight of 670 pounds has been assumed for the above 
satellite. The center of mass is assumed to be at the 
center of the rectangular prism. An attitude reference 
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The total injection sequence from the launch configura- 
tion to mission mode is now briefly presented. The launch 
into a near polar orbit, is assumed to be retrograde out of 
Cape Kennedy, and orbit is achieved by direct ascent on 
the descending node approximately eighteen minutes after 
lift-off. Within the next two minutes, the satellite is de- 
tached from the third stage by a force derived from mul- 
tiple separation springs. Immediately after separation, 
the total satellite is despun, by release of its own yo-yo 
weights, from the orbit-injection third stage spin rate of 
ninety revolutions per minute to a nominal 3.2 revolutions 
per minute which provides the design value of momentum 
for the satellite attitude control system. After this rapid 
spin-rate reduction, the control system flywheel is ener- 
gized to a fixed 75 revolutions per minute, which reduces 
the satellite spin rate to 1.6 revolutions per minute by a 
transfer of momentum. In this mode, the satellite is com- 
pletely spin-stabilized about the flywheel (launch-direction) 
axis, and a magnetic torquing procedure can be employed 
for turning the momentum vector (and spin axis) from the 
injection attitude and aligning it with the orbit normal. 

When the spin axis has been moved to within ten degrees 
of the orbit normal, the pitch-axis control can be com- 
manded to achieve local vertical orientation of the satellite, 
transferring most of the total momentum of the satellite 
to the flywheel. The satellite is then rotating at only one 
revolution per orbit and the flywheel at approximately 150 
revolution per minute, and the solar panels can be erected 
to complete the acquisition sequence. 


Satellite Orbit and Injection Errors 


It is known that an earth satellite orbit which is inclined 
to the equator will undergo a gyroscopic precession due to 
the earth's equatorial bulge. This precession is easily ob- 
tained from Equation 2-d repeated from the previous sec- 
tion as* 
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♦The total effect on an earth satellite, of the gravitational 
potential for an oblate body having axial symmetry, has 
been analyzed in a far more detailed and rigorous man- 
ner in Chapter XVII of the Reference 1 . 


Figure 1. 670 pound sun-synchronous 
meteorological satellite. 
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Figure 2. Curve defining sun-synchronism. 


The perturbing acceleration due to the equatorial bulge 
has been treated extensively in the literature (see for 
example Reference 4). The component, W , normal to the 
satellite orbit is 

W = - 3J 20 sin I cos I sin(uj + 6) (8) 


Substitution into Equation 7 provides 
, 2 
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For a circular orbit, the average precession over one 
orbit becomes 
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For the satellite to be sun-synchronous its orbit must 
precess at the average rate of 0.985 degrees per day, the 
rate at which the earth revolves about the sun. Figure 2 
is the result of setting Equation 10 equal to this rate. Any 
combination of altitude and inclination falling on this curve 
defines a sun-synchronous orbit. A relatively low altitude 
of 750 nautical miles and the corresponding inclination of 
101.4 degrees have been chosen for the reference orbit in 
for this study. 

Since the precession rate is a function of altitude and 
inclination, an injection error in achieving the nominal 
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Figure 4. Precession error vs. altitude error. 


parameters will cause a variation from the desired rate of 
precession. Figures 3 and 4 show the effect on precession 
of altitude and inclination variations, respectively (around 
the nominal values of 750 nautical miles and 101.4 degrees) 
The precession error shown is the difference between the 
actual precession and the desired rate. The total error in 
the desired rate of precession is the sum of the effects of 
the altitude error and inclination error. For both altitude 
and inclination errors, any resultant positive "drift" in 
precession is in the eastward direction. 


Shown as Table 1 are estimates of the one sigma (lcr ) 
booster errors in achieving the nominal parameters.* 
These are also indicated as points on the plots of Figures 
3 and 4. Note that inclination error is the source of nearly 
all of the root-sum-square error of 0.029 degrees per day 
precession rate error. 

Orbit Correction 


Because of the injection errors described in the pre- 
ceding section, a means of orbit correction must be pro- 
vided for sun-synchronous satellites. The precision of 
on-board equipment will in many cases dictate the method 
of orbit correction which must be employed. For purposes 
of this paper, however, it is assumed that any injection 
within the limits of error herein provided is satisfactory 
for on-board equipment operation. Therefore, orbit cor- 
rection is necessary only to maintain sun-synchronism. 


Table 1 

Estimated One Sigma ( lcr) Booster Errors in 
Achieving Nominal Injection Parameters. 


One Sigma 

Precession Error 

Error 

(Degree/Day) 

Inclination — 0.33 Degrees 

0.028 

Altitude -10.8 N. MI. 

0.008 


RSS - 0.029 


Figure 3. Precession error vs. inclination error. 


*W. R. Schindler to Herman I.agow, U.S. Government 
Memorandum, dated October 21, 19G6, "Tiros M Sys- 
tems Review." 
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The method of correcting to the sun-synchronous curve 
by changing inclination can be analyzed by Equation 6, re- 
peated here as 

SI = - (F/m) St. (13) 

77 V 

c 

This equation can be rearranged to 

I T = FSt = g) mV c Si (14) 

The maximum injection error of 0.029 degrees per day in 
precession is equivalent to an inclination error of 0.34 
degrees. Therefore, 4470 pound-seconds of total impulse 
is required for orbit correction by inclination adjustment. 

Analysis of the method of correcting to the sun- 
synchronous curve by changing altitude must be approached 
differently. First, Equation 3 can be rearranged and 
written as 


Figure 5. Injection with respect to nominal orbit. 


This being the case the method of correction chosen 
should be that which minimizes propulsion system re- 
quirements. 


Figure 5 shows a typical injection, altitude and inclina- 
tion, of a sun-synchronous satellite with respect to the 
nominal condition. One method of maintaining sun- 
synchronism is to continuously thrust in a manner which 
would cause a precession equal in magnitude but opposite 
in direction to that (a maximum of ±0.029 degrees per 
day) resulting from the injection errors. Another method 
is to correct either or both the altitude and inclination to 
the curve defining sun-synchronism. In this paper cor- 
rection of altitude and inclination will be considered sep- 
arately and compared. 


Equation 5 repeated here as 

sin I S (—] = — ^ (ID 
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expresses the necessary relationship for continuous thrust 
and precession. After rearrangement the total impulse 
is: 
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which is correct to first order. Substituting Equations 15 
and 17 into Equation 16, retaining first order variations 
only and rearranging yields 
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The indicated calculation using previously suggested 
nominal values gives a total impulse requirement of 2050 
pound -seconds. 


Table 2 summarizes the total impulse requirements for 
correcting the precession rate to the sun-synchronous rate 
As previously stated, the method of continuous thrust and 
precession requires an excessive total impulse. Of the 
two remaining methods, correcting altitude is less expen- 
sive, less than one-half, than correcting the inclination. 


Table 2 


In order that a comparison can be made with the other 
methods a requirement for one year's maintenance (St = 
1 year) is established. Substitution of the maximum pre- 
cession error of 0.029 degrees per day for S (dfi/dt ) 
and the nominal values for the remaining terms yields a 
total impulse requirement of 139,000 pound-seconds. As 
will be seen in the following paragraphs, this total impulse 
is excessive when compared to that required by the other 
methods of correction. 


Comparison of Methods of Correcting to 
Sun-Synchronous Precession Rate. 


Method of 

Total Impulse 

Correction 

(Pound-Seconds) 

Altitude Adjustment 

2,050 

Inclination Adjustment 

4,470 

Constant Precession (One Year) 

139,000 





Therefore, the method of correcting altitude is chosen 
for the remaining analysis of the sun-synchronous satel- 
lite. However, an effect on the orbit which may be im- 
portant should be noted. It is remembered that of the al- 
titude and inclination injection errors, the inclination 
error is the major contributor to the precession rate 
error (see Table 1). Thus, the possibility of having a 
final orbit as much as 35 nautical miles above or below 
the original nominal value, Figure 5, exists because an 
inclination error is being corrected by changing altitude. 
This is equivalent to changing the orbital period by ap- 
proximately 1.5 minutes. For present purposes these 
effects are assumed not to be detrimental. 

An especially attractive feature of performing the orbit 
correction by changing altitude is that it requires tangen- 
tial thrust. As for a number of meteorological missions, 
the satellite is required to be stabilized about the roll and 
yaw axis, and motion about the pitch axis is required to 
be controlled to provide an earth oriented platform. The 
two thrusters, one for raising and one for lowering the 
orbit, can therefore be rigidly mounted to the satellite 
such that they are always aligned tangent to the orbit 
when the satellite is in its desired attitude. The thrust- 
ers are mounted 20 inches from the center of mass. 
Because of imprecise knowledge of center of mass loca- 
tion and other factors, misalignment errors in mounting 
the thrusters to the satellite body exist. In this paper the 
thrusters are assumed to be aligned within a one degree 
cone angle with respect to the center of mass. 

Attitude Control Limitation 

The attitude control system assumed, References 2 and 
3, has been used for a number of meteorological satellite 
missions. Reaction between the earth's magnetic field 
and a magnetic field generated within the satellite pro- 
vides control over yaw and roll components of attitude 
and the total system angular momentum. A sampled data 
servo system with pitch error input (derived from atti- 
tude relationship between the satellite and the earth's 
horizon) provides the required control of the satellite 
pitch motion by means of momentum interchange between 
the wheel and satellite. 

By its nature, a control system of this type has limited 
capability. It has been estimated that the maximum al- 
lowable deviation of the spin vector from its desired posi- 
tion is ten degrees per orbit.* Anything in excess of this 
amount could be beyond control limits. Therefore, the 
maximum angular momentum disturbance causing spin 
vector motion is 35 inch-pound-seconds corresponding to 
the satellite and flywheel angular momentum of 200 inch- 
pound-seconds. The spin momentum removal capability 
has been estimated to be 2.5 percent or five inch-pound- 
seconds per orbit. 

Thruster Selection and Operation 

The requirements of the sun-synchronous satellite 
mission described in this paper are such that either a pair 
of ion engines or a pair of resistojets are capable of per- 
forming the orbit adjustment. For purposes of compari- 
son a 300 micropound ion engine and a five millipound 

* These estimated control system limitations are based on 
data taken from Reference 3. 


resistojet were selected. Both of these systems have 
weight and power requirements of approximately 25 pounds 
and 70 watts respectively. Two other important consider- 
ations are (1) time required to make the correction and 
(2) operation of the thruster within the control system 
limits. The latter must of necessity be considered first. 

In an earlier section the thrusters were located 20 
inches from the center of mass at a possible maximum 
misalignment angle of one degree. Anytime the thruster 
is on, a disturbance torque of 

T d = 20F tan 1 ° = 0.35F inch-pounds 

where F is thrust level, is applied to the satellite. The 
angular momentum introduced over an interval of time 
h t is therefore 

h d = T d St = 0.35 F St inch-pound-seconds. 

Earlier it was established that h , should not exceed five 

d 

inch-pound-seconds per orbit. Substituting this along with 
the five millipound thrust level for the resistojet into the 
above equation shows that the maximum disturbance level 
is exceeded in 47.6 minutes or less than one half the 
orbital period of 113.5 minutes. The resistojet can there- 
fore be operated for only 47.6 minutes of each orbit. How- 
ever, a similar calculation shows that the ion engine can 
be operated continuously. 

Having calculated the operational time per orbit, the 
total correction time can be calculated from 

I ( = F St = 2050 pound-seconds 

which is the total impulse required for the tangential cor- 
rection. The calculation shows that 11.3 days is required 
by the resistojet and 79.1 days by the ion engine. Table 3 
summarizes the system characteristics and operation for 
the resistojet and ion engine. A more detailed description 
of a typical propulsion system will be presented for the 
earth-synchronous satellite. 

Earth-Synchronous Satellite 
Satellite Description 

A 1400 pound communications satellite, Figure 6, was 
chosen as typical of a class of earth-synchronous satel- 
lites to be flown within the next decade. The satellite con- 
sists of two equipment modules separated 15 feet by a 


Table 3 

Summary of Characteristics for Orbit 
Correction Propulsion Systems. 



Ion Engine 

Resistojet 

Thrust (Millipounds) 

0.3 

5.0 

Weight (Pounds) 

25 

25 

Power (Watts) 

70 

70 

Time Required (Days) 

79.1 

11.3 

Operation Per Orbit 

Continuous 

47.6 Minutes 
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Launch Profile 


z 



Figure 6. 1400 pound earth- synchronous 
communications satellite. 


six-member truss. All satellite subsystems and com- 
ponents are housed in these modules. One equipment 
module is earth viewing; the other views space and is 
referred to as the "aft equipment module." A thirty-foot- 
diameter parabolic antenna reflector is mounted to the 
aft equipment module. Solar cell paddles for the solar 
conversion power supply extend beyond the rim of the re- 
flector. They are mounted on two extended and reinforced 
antenna deployment trusses. One panel is rotated 90 de- 
grees with respect to the other to provide greater aver- 
age solar radiation incidence per orbit. 

A reasonable center of mass location for this config- 
uration is six feet from the space viewing face of the aft 
equipment module on the axis of symmetry of the para- 
bolic dish. For later reference a coordinate system, 
rigidly fixed to the satellite and with origin at the center 
of mass is defined as follows: 

X-axis or yaw axis — coincides with the axis of sym- 
metry of the parabolic dish. 

Y-axis or roll axis — is normal to the plane containing 
the solar panel booms and yaw axis. 

Z-axis or pitch axis — lies in the plane containing the 
solar panel booms and yaw axis. 

Representative moments of inertia for this configuration 
are: 

I x = 1850 slug -ft 2 , 

I = 2690 slug-ft 2 , and 

y 

I 2 = 1620 slug-ft 2 . 

The desired attitude of the satellite is to have the yaw 
axis of the satellite pointing to some point on the earth's 
surface and to have the roll axis near the satellite 
velocity vector. 


For purposes of this paper it is assumed that the satel- 
lite has been successfully injected into a near synchronous 
orbit (altitude equals 19,323 nautical miles) with residual 
injection errors as given at the end of this section. Those 
aspects of the launch profile leading up to this injection 
are now briefly reviewed. 

The launch profile is assumed to be similar to that of 
existing operational satellites (e.g. Syncom, Early Bird, 
and ATS-B); i.e., the satellite/launch vehicle combination 
injects into a low (90 to 100 nautical miles) circular in- 
clined orbit and then at some equator crossing the satel- 
lite is injected into a Hohmann transfer ellipse having a 
perigee altitude equal to that of the low circular orbit and 
an apogee altitude equal to the synchronous altitude. At 
some apogee passage of the transfer orbit a propulsion 
system attached to the satellite is ignited, injecting the 
satellite into a near circular equatorial orbit, i.e., the 
thruster simultaneously removes the eccentricity and in- 
clination of the transfer orbit. 

Because of non-perfect systems performance, injec- 
tion errors are to be expected. For present purposes, 
estimated three sigma (3o-) dispersions for launch vehicle 
and spacecraft — induced errors were assumed based on 
previous missions and mission studies similar to the cur- 
rent one. The five sources of error considered were 
(1) induced 3 a errors for the transfer to the synchronous 
orbit, (2) total impulse dispersion of the apogee kick motor, 
(3) attitude pointing error, (4) net thrust reduction result- 
ing from vehicle coning during the firing of the apogee 
motor, and (5) a timing error in igniting the apogee motor. 
It should be noted that spin stabilization has been assumed 
for the transfer orbit. Table 4 contains the estimated root- 
sum-square errors from the above sources. From these 
errors, the in-plane boundaries within which the actual 
orbit should fall can be calculated. The results are shown 
as Figure 7. 

Since it is desirable to remove eccentricity and inclina- 
tion with as high a thrust as possible to keep pulse time to 
a minimum and since spin stabilization has been assumed 
for the transfer orbit (this implies a requirement for pro- 
cessing a large angular momentum) a relatively high thrust 
is consistent with removingthe above errors. A five pound 
hydrazine system aligned to thrust along the spin axis 
should satisfactorily perform these functions. It is esti- 
mated that sixteen hours after the injection into the syn- 
chronous orbit, sufficient data will have been obtained to 
provide a preliminary estimate of where within the bound- 
aries of Figure 7 the actual orbit lies, and that a final orbit 
determination can be made within 28 hours. Based upon 

Table 4 

Estimated Root-Sum-Square Injection Errors 
into a Synchronous Orbit. 

Altitude - ± 100 nautical miles 

Velocity - ± 90 feet per second 

Inclination - ±0.5 degrees 

Eccentricity - 0.022 
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Figure 7. In-plane orbital error boundaries 
due to injection errors. 

this information and the further assumption that the orbit 
injection point must be one of the nodes, if there is any 
inclination, the sequence of precession and orbit correc- 
tion can be predicted. Depending on whether the actual 
orbit is high or low as depicted in Figure 7, the sequence 
is as shown in Table 5. It is noted that the satellite has 
been assumed to remain in the attitude required for syn- 
chronous orbit injection until the first maneuvers indi- 
cated in Table 5 are made. It is further noted that the 
inclination removal is made only after sufficient time for 
an accurate orbit determination as no inclination error 
should remain after this correction. However, the pre- 
liminary orbit determination has been assumed sufficiently 
accurate for eccentricity correction of the low-orbit con- 
dition as additional in-plane orbit adjust capability is 
being provided by the low thrust propulsion system. Since 
all in-plane velocity increments have been applied at or 
near the synchronous orbit injection point (see Table 5), 
the altitude error of ±100 nautical miles at injection still 
remains. So that the proper orbit period is achieved, any 
error at injection should be offset by an equal and oppo- 
site error 180 degrees from injection, i.e., if the orbit is 
100 nautical miles high at injection it should be made 100 
nautical miles low 180 degrees from injection. There- 
fore, the orbit correction remaining for the low thrust 
system lies somewhere between two circular orbits of 


Figure 8. In-plane orbital error boundaries 
of errors to be removed by low thrust. 

radii 

r = r ±100 nautical miles 

S 

as depicted by Figure 8. 

Low Thrust Propulsion System Requirements 

A brief description of the launch sequence has just been 
given. It is now assumed that despin, down to one degree 
per second body rates, has been accomplished by yo-yo's 
and that the parabolic antenna and solar paddles have been 
deployed. This leaves the satellite in its final flight con- 
figuration at low body rates. Functions left still to be ac- 
complished are (1) initial attitude acquisition, (2) vernier 
orbit adjustment, (3) operational attitude control, includ- 
ing fine pointing and slewing maneuvers, and (4) station 
keeping. 

The initial attitude acquisition might consist of estab- 
lishing a lock of the satellite roll axis on the sun line, 
i.e., control pitch and yaw position and roll rate. The 
local vertical can then be acquired by rolling around the 
sun line. Finally, acquisition of Polaris would establish 
the yaw attitude. The performance of these maneuvers 
requires sufficient thrust for completion in reasonable 
time without exceeding limiting angular accelerations of 
the antenna in the deployed condition. These angular 


Table 5 

Sequence of Removing Eccentricity, "e", and inclination, "I", 
with the Hydrazine Propulsion System. 


Low Orbit 

High Orbit 

Injection point is 
descending node 

Injection point is 
ascending node 

Precess 20.7° 

Precess 69.3° 

Precess 110.7° 

♦Remove "e" at 24 hrs. 

♦Remove "1” at 36 hrs. 

♦Remove "I" at 36 hrs. 

Precess 90° 

Precess 90° 

Precess 90° 

♦Remove "I" at 36 hrs. 

♦Remove "e" at 48 hrs. 

♦Remove "e" at 48 hrs. 


♦Time is measured from injection into the synchronous orbit. 
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accelerations, assumed to be ten degrees per minute 
squared for pitch and roll and five degrees per minute 
squared for yaw, result in the following limiting torque 
levels. 

Roll torque — 0.13 ft-lb 
Pitch torque — 0.08 ft-lb 
Yaw torque — 0.04 ft-lb 

These torques will be used later for sizing the thrusters. 
The total angular momentum, and therefore fuel, required 
of the thrusters for acquisition is small compared to that 
required for operational attitude control and will be in- 
cluded in the latter. 

Vernier orbit adjustment is the second of the functions 
required of the low thrust system. At the completion of 
the orbit corrections described as part of the launch pro- 
file, no inclination error should remain. However, the 
in-plane orbital errors bounded by 

r = r ± 100 nautical miles 

S 



Figure 9. Self optimizing adaptive control for 
low disturbance torques and minimum control. 



as shown by Figure 8 still remain. A total velocity in- 
crement of 22.0 feet per second is required to circularize 
the worst case orbit which is 100 nautical miles low at 
perigee and 100 nautical miles high at apogee. It is noted 
that this velocity increment is based on optimum timing, 
meaning that thrust should be applied at or near (say 
within ± 1.5 hours) apogee or perigee. In order to main- 
tain temperature, and therefore specific impulse, the 
thruster is operated in a duty cycle mode. During each 
pulse heat is lost through the propellant flow, reducing 
temperature. A ten percent duty cycle, which is assumed 
for present purposes, provides a sufficient coast period 
for the original temperature to be regained. As will be 
seen later the tangential component of thrust to be applied 
to the 1400 pound satellite for three hours twice a day is 
0.0174 pounds. The total correction time is then 26 days 
for the worst case. 

Operational attitude control must be maintained for 
both fine pointing and slewing modes of operation. A 
requirement for slewing from horizon to horizon (a total 
of 17.5 degrees) and being restabilized within thirty min- 
utes seems reasonable. This places no additional hard- 
ship beyond that of the initial acquisition on the control 
system. The angular momentum required for each pitch 
slew is 0.94 foot-pound-seconds and for each roll slew 
is 1.58 foot -pound-seconds. The coupling torque on the 
yaw axis due to a pitch/ roll slew (one degree/minute) is 
about 10~ 4 foot-pounds which represents about 0.1 foot- 
pound-seconds per slew. Assuming one pitch and one 
roll slew per day, the total angular momentum required 
for two years is 1910 foot-pound-seconds. 

A fine pointing accuracy of 0.1 degrees is assumed to 
be maintained by a self -optimizing and adaptive control 
logic technique based on the use of a small special pur- 
pose computer to monitor and update the system gain to 
improve performance. Ideally, during periods of a given 
mission when external bias torques are very small, it is 
desirable to make the minimum control torque impulse 
as small as possible, thereby minimizing the frequency 
of attitude deadband crossings (and therefore minimizing 


Figure 10. Self optimizing adaptive control for high 
disturbance torques. 

fuel consumption). However, during those periods of the 
mission when a relatively high bias torque exists, a para- 
bolic attitude trajectory will occur in the deadband after 
each unidirectal torque impulse. These pulses should be 
just large enough to cause the trajectory to cross the 
deadband but not touch the opposite threshold, thus de- 
creasing the number of thruster firings (increasing 
thruster reliability) and increasing specific impulse be- 
cause a higher temperature can be maintained. These 
limit cycle performance characteristics are illustrated 
as Figures 9 and 10. 

The angular momentum required to compensate for the 
solar pressure disturbance is the integral of the absolute 
value of disturbance torque. The solar pressure torque 
profile has been assumed to be sinusoidal with a maximum 
value of 10~ 4 foot-pounds for both the pitch and roll axes. 
Therefore, 

Angular Momentum 
Period 

4 x IQ- 4 

CO 

where o, is the frequency of the disturbance. The period 
for the pitch disturbance is the orbital period and is one 
year for the roll disturbance. The total angular momen- 
tum required for two years fine pointing is therefore 8010 
foot-pound-seconds. Therefore for slewing and fine point- 
a total of 9920 foot-pound-seconds angular momentum 
must be provided. For purposes of later sizing the pro- 
pulsion system a requirement for 12,000 foot-pound- 
seconds is assumed. This should cover initial acquisition 
and provide for additional operations such as, for example, 
tracking a low orbit satellite. 

Station Keeping is the fourth and final task to be per- 
formed by the low thrust propulsion system. From this 


2 x 10 


J -P/ 2 
n 


sin &Jt dt 
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standpoint, 53° West Longitude, which can be reached on 
the second pass of apogee in the transfer orbit, is the 
worst-case condition as the J 22 gravitational perturba- 
tion is maximum at this point. This station could theoret- 
ically be maintained by applying a constant, very small 
thrust against the perturbation which would produce a 
total velocity increment of 5.36 feet per second per year. 

Thruster Configuration and Sizing 

The satellite configuration is such that the thrusters 
cannot be located at the center of mass. In addition a 
constraint is assumed requiring maintenance of attitude 
control even with one thruster inoperative. Under these 
conditions a thermal storage resistance multijet system, 
to be described in detail in the following section, is cap- 
able of performing the required functions. 

The thruster configuration is shown as Figure 11. All 
of the thruster jets are located in the same plane, per- 
pendicular to the yaw axis and 4.6 feet above the satellite 
center of mass. Each thruster is located 3.6 feet from 
the yaw axis. This configuration satisfies all of the above 
requirements providing thrust is applied simultaneously 
through two nozzles, e.g. TM21 and TM23, to give a re- 
sultant vector through the center of mass for the orbit 
adjustment and station keeping functions. This technique 
can be successfully applied because the radial component 
of thrust has negligible effect on the orbital elements 
when applied in the vicinity of apogee and/or perigee, see 
Equations 2-a and 2-b, and therefore only provides the 


TM33 



Figure 11. Thruster configuration. 


component of thrust necessary to balance the torque ap- 
plied to the satellite by the station keeping component. 

The requirements of all functions to be performed by 
the thrusters must be taken into account when sizing each 
jet. It is remembered that there must be sufficient thrust 
to acquire the sun and earth in reasonable time during the 
acquisition phase. Also, a small impulse for attitude con- 
trol during the local vertical fine pointing limit cycle is 
required. Investigations showed that the minimum pulse 
(limited by jet "on" time) could be easily met. For this 
reason it was decided to use as high a thrust level as pos- 
sible. The maximum jet size, therefore, is limited by the 
maximum specified acceleration the parabolic antenna can 
sustain in the deployed condition. The resulting torque 
levels, given earlier as 0.13, 0.08, and 0.04 foot-pounds 
for roll, pitch, and yaw respectively, lead to the following 
thrust levels: 


TM11 = 0.0361 pounds 
TM12 = 0.0056 pounds 
TM13 = 0.0283 pounds 
TM21 = 0.0222 pounds 
TM22 = 0.0056 pounds 
TM23 = 0.0174 pounds 

Symmetry dictates that modules TM1 and TM3 and modules 
TM2 and TM4 have the same thrust levels. Table 6 gives 
the various propulsion functions of each thruster, includ- 
ing attitude control with one module imperative. 

Thermal Storage Resistance Jet 


As discussed in the previous section, a thermal stor- 
age resistance multijet system has been chosen to pro- 
vide for the low thrust station keeping and attitude con- 
trol requirements. Ammonia is the selected propellant 
for this system. Figure 12 shows the essential features 
of a typical resistance jet. The overall length is about 
five inches with approximately a three inch diameter. The 
thruster module weighs less than 1.5 pounds. The thruster 
assembly consists of a body with four nozzles in its outer 
shell, a heater subassembly consisting of resistance wire 
wound on an electrical insulator, and heat shielding fabri- 
cated from layers of thin metallic foil separated by mini- 
mum contact support wires. 


The purpose of the resistance jet is to increase the 
specific impulse of the gas by the addition of heat. As the 
propellant passes through the thruster body, it is heated 
up to operating temperature by its contact with the hot flow 
passage. This hot surface receives its heat from the 
heater unit which is installed in the center of the thruster 
body. 

The resistance jet can be operated at temperatures up 
to 2000°F. This operational flexibility presents the option 
of trading off specific impulse and power input. Figure 13 
shows an estimate of the resistance jet specific impulse 
to power input for short pulses. This curve is based on 
test data from existing thrusters together with estimates 
of design improvements now being incorporated. 
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SPECIFIC IMPULSE (SECONDS) 


Table 6 


Resistojet Propulsion System Utilization. 


Thruster Jet 
Number 

TM11 

TM12 

TM13 

TM21 

TM22 

TM23 

TM31 

TM32 

TM33 

TM41 

TM42 

TM43 

Jet Thrust Level 
(Millipounds) 

36.1 

5.6 

28.3 

22.2 


17.4 

36.1 

5.6 

28.3 

22.2 

5.6 

17.4 

Station Keeping 




* 


* 







+ Roll Torque 



** 




* 






- Roll Torque 

* 








** 




+ Pitch Torque 




* 








** 

- Pitch Torque 






** 




* 



Yaw Torque 
(Couple) 


** 






- 

** 



* 



♦Signifies the jets to be used with all four modules operative. 
♦♦Signifies the jets to be used if the * jets had failed. 


THERMAL 

SHIELDING 



Jet nozzle throat sizes and contours are designed for 
a particular operating condition. The thrust level of the 
jets is dependent on the nozzle design and both the final 
temperature and pressure of the gas. Existing designs 
range from about five to 60 millipounds thrust with supply 
pressures ranging from less than one up to approximately 
five atmospheres. The gas temperature will depend on 
the power input to the heater. 

The thrust level of the resistance jets takes a finite 
amount of time to build up to the full thrust level. Tests 
run to date show that it takes approximately 0.1 second to 
reach approximately 95 percent of full thrust. It also re- 
quires a finite time for the thrust level to decay from full 
thrust down to zero. The build up and decay time will be 
dependent on the solenoid reaction time as well as the 
volume of propellant between the solenoid and the thruster 
nozzle. 


Figure 12. Typical thermal storage 
resistance multi-jet. 



0 4 8 12 16 20 24 28 32 36 

POWER (WATTS) 

Figure 13. Estimated short pulse specific impulse 
vs. power input for a thermal storage resistance jet. 


Propellant Required 

In a long-life mission with fine pointing requirements 
one of the primary concerns is fuel consumption. It is re- 
membered that 22 feet per second was required for vernier 
orbit adjustment and 5.36 feet per second per year or ap- 
proximately 11 feet per second for a two year mission for 
station keeping. Both of these numbers must be multiplied 
by the factor (3.6 + 4.6)/3.6 to account for the radial com- 
ponent of thrust needed to balance the body torques. The 
total translational velocity requirement is therefore 75 
feet per second. For the 1400 pound satellite this means 
a momentum requirement of 3260 pound-seconds and, 
assuming a fuel specific impulse of 200 seconds, requires 
16.3 pounds of fuel. The total angular momentum required 
of the propulsion system was estimated earlier to be 12,000 
pound -feet-seconds. If a four foot moment arm is assumed 
(actually 3.6 feet part of the time and 4.6 feet the remainder 
der) and again assuming a specific impulse of 200 seconds, 
15.0 pounds of fuel is required for initial acquisition and 
operational attitude control. A total of 57.0 pounds of fuel 
is being specified to provide a contingency of 25.7 pounds. 
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Figure 14. Ammonia feed system schematic diagram. 


Ammonia Feed System 

In order to obtain a high degree of accuracy in predic- 
ting thrust levels of the resistance jets, a constant pres- 
sure zero gravity storage and feed system is used to pro- 
vide the ammonia propellant for the thrusters. The 
system chosen for this application is shown schematically 
in P'igure 14. The storage tank shown is approximately 
17.7 inches in diameter with a 7.6 inch regulator assembly 
attached to one end. Fully loaded, the tank can hold 57 
pounds of ammonia. The ammonia in the storage tank is 
maintained at the saturation pressure (up to a maximum 
of about 200 psia) corresponding to the local environ- 
mental temperature. Including the fuel, the maximum 
weight for this subsystem would be 87 pounds. 

The regulator assembly attached to the tank provides 
vapor flow at a regulated pressure within one percent of 
a preselected value (between one and five atmospheres). 
This regulator satisfies the zero gravity requirement: 
vapor flow output with either liquid or vapor flow from 
the .storage tank. 

Propellant flow into the resistance jet is accomplished 
by applying a signal to the solenoid valve adjacent to the 
f m uster. As the gas leaves the plenum tank the pressure 
in both the plenum and preplenum drops. The pressure 
drop in the preplenum closes the pressure switch which 
opens the valve providing a pulse of liquid or vapor am- 
monia from the storage tank into the preplenum. Any 
necessary vaporization takes place in the preplenum. The 
vapor then flows through an orifice which eliminates 
deadband and overshoot conditions in maintaining the 
plenum pressure within one percent of the preselected 
value. 

For system reliability, the regulator assembly is 
equipped with a filter, a rupture disk and relief valve as- 
sembly, and a redundant pressure switch and solenoid 
valve. The storage tank also contains a repture disk and 

relief valve assembly. 

The total weight of this system is 103 pounds. This 
includes 57 pounds of propellant, ten pounds for four 
power conditioners, six pounds for the four thrusters, 
ai.ci : 0 pounds for the storage and feed system. 

Conclusions 

Analyses of applications of electric propulsion to sun- 
v't/vhronous and earth-synchronous satellite missions 


have been performed. These analyses were made using a 
perturbation form of the equations of motion first used by 
Gauss . 

For the sun-synchronous satellite, it was necessary to 
provide orbit correction of injection errors causing an 
orbit precession error of 0.029 degrees per day. It was 
found that: 

• Tangential thrusting (raising or lowering altitude) 
was a more efficient method of correction than 
changing inclination. Providing precession equal to 
but opposite in direction to the error was found to 
require an excessive total impulse. 

• Either a 300 micropound ion engine or a five milli- 
pound resistojet can be used for removing injection 
errors. Both systems weigh approximately 25 
pounds and require approximately 70 watts of power. 
The resistojet required 11.3 days for the correction, 
more than half of which is spent coasting due to con- 
trol system limitations. The ion engine can be 
thrusted continuously but still required 79.1 days be- 
cause of the low thrust magnitude. 

It was found that a thermal storage resistance multijet 
system could provide for the initial attitude acquisition, 
vernier orbit adjustment, operational attitude control, and 
station keeping requirements. With respect to these four 
functions it was found that: 

• Limiting torques of 0.13, 0.08, and 0.04 foot-pounds 
respectively in roll, pitch and yaw during initial 
acquisition determined the thruster nozzle sizes. 
Thrust requirements varied between 5.6 and 36. 1 
millipounds. 

• A 22 foot per second vernier velocity increment was 
required to circularize the orbit. Application of 
17.4 millipounds of tangential thrust for three hours 
twice a day for 26 days was required. 

• Operation attitude control, both for fine pointing and 
for one slew maneuver per day for two years, re- 
quired 12,000 foot-pound-seconds total angular 
momentum. 

• Station keeping required application of 5.36 foot per 
second per year through the 17.4 millipound nozzles. 

Ammonia was used lor the propellant. The propulsion 
system was found to weigh 103 pounds of which 57 pounds 
was ammonia. 
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